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Resume. We study experimentally the paths of an assembly of cracks growing in 

interaction in a heterogeneous two-dimensional elastic brittle material submitted to 

uniaxial stress. For a given initial crack assembly geometry, we observe two types 

of crack path. The first one corresponds to a repulsion followed by an attraction on 

^ I one end of the crack and a tip to tip attraction on the other end. The second one 

S ' corresponds to a pure attraction. Only one of the crack path type is observed in a 

1-^ , given sample. Thus, selection between the two types appears as a statistical collective 

5Ih ' process. 
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1. Introduction 

Understanding tlie rupture mechanisms of solids has become an important goal of 
fracture physics in order to improve the strength of structures and avoid catastrophic 
failures. Characterization of rupture properties very often involve the growth of a 
dominant crack. For instance, there is an extensive literature, both experimental and 
theoretical, discussing the slow growth dynamics of a single crack along a straight path 
in brittle [BElEllllElElITj or visco-plastic materials [HlElIinillllllllllSlITllIS]. 
However, in many practical situations, the crack path is not straight. For one thing, 
the crack path can be slightly destabilized and develop some roughness due to 
dynamical instabilities [16] or material heterogeneities. This kind of path instability 
has motivated a large amount of experimental works analyzing the roughness of cracks 
[T71 fT8l fT9l [20l [211 [221 [23] as well as several models describing roughening mechanisms 
[M] [25] [26] and non-trivial effects of heterogeneities on the rate of crack growth 
[271 [28l [29] . Much larger deviations of the crack path from a straight line can be 
observed during the growth of an array of interacting cracks. Comparatively to the case 
of roughness, this rather complex situation has been studied essentially theoretically 
[ani [an [321 [331 [31 [351 [36] and very little experimentally [371 [381 [39] despite its practical 
importance especially for heterogeneous materials where multiple cracks are likely to 
form. Understanding the growth of interacting cracks is also a relevant issue for fault 
dynamics [IQ] as well as crack pattern formation during drying processes [HI [121 [l3l [M] . 

There are several levels of complexity that can play a role in the growth of 
interacting cracks. First, depending on the geometry of the crack array and the loading, 
the stress field around a crack will be amplified or shielded because of the existence 
of other cracks. Then, in 2d, the effective contribution of the loading on each crack 
will usually be a mixture of mode I and mode II |j|. It has been argued that crack 
deviations occur when the shear stress on the crack lips (mode II) is non-zero and 
that the crack will grow so as to minimize the mode II contribution and maximize the 
mode I contribution ^5l [46l [47]. A noticeable observation related to this property is 
the preferential merging of cracks joining each other forming a right angle, commonly 
observed for crack patterns in drying experiments [121 [131 HI]. This occurs naturally 
because the principal stress along a crack lip is parallel to the crack direction so that 
a crack approaching the lip with a right angle is propagating mainly in mode I. On 
the other hand, when two collinear mode I cracks are growing towards each other, they 
do not merge tip to tip, but instead repel each other [37]. The origin of this effect has 
been discussed by several authors [30l [32] and numerical simulations have been able to 
reproduce, at least qualitatively, experimental observations [3l] . However, there are not 
many experimental data to compare with these theoretical predictions, and furthermore, 
there is very little knowledge about the effect of material heterogeneities on the crack 
path selection. 

X The mode I corresponds to a tensile stress normal to the plane of the crack as the mode II corresponds 
to a shear stress acting parallel to the plane of the crack and perpendicular to the crack front. 
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In this paper, we study experimentally the trajectories of interacting cracks in an 
almost two dimensional brittle and heterogeneous material. The array of cracks has 
been initiated prior to the application of the external stress and the heterogeneity 
of the material (paper sheets) allows us to test the stability of crack paths to small 
perturbations. We uncover that for a given geometry of the crack array, a crack can 
follow statistically two stable trajectories : an attractive one and a repulsive one with 
respect to the neighboring cracks. The main result of this investigation is the analysis 
of the geometrical conditions for which cracks are attracted towards another and when 
they are repelled. The paper is organized as follows. In the next section, we describe 
the experimental apparatus. In section 3, the extraction procedure of the post-mortem 
shapes of the crack profiles is described. In sections 4, 5 and 6, we analyze the different 
types of crack paths and their statistics as a function of the initial crack pattern 
geometry. In section 7, we finally discuss the results and conclude. 

2. The experiments 

To study the growth of several interacting cracks, we have loaded bi-dimensional 
brittle samples made of fax paper sheets (Alrey®) with a tensile machine. The samples 
were previously prepared cutting several straight cracks in the paper sheets. In order 
to observe clearly the interaction between cracks during their growth, we had to find 
a geometry for which all crack tips are equivalent. Actually, we wanted to prevent the 
isolated growth at a particular crack tip that will inhibit the growth of the other cracks. 
One possible geometry that follows this condition consists in an array of cracks as 
presented in figure [H The pattern formed by the two lines of cracks offers the advantage 
to have some translational invariance in the cracks direction. The stress is applied 
perpendicularly to the cracks direction and uniformly on the sample borders. Therefore, 
stress intensification is theoretically equivalent at each crack tip. In this geometry, there 
are three adjustable parameters : the crack length, the vertical spacing of cracks on 
a line and the horizontal spacing between the two lines. In this article, the first two 
parameters are fixed to 1cm and only the distance d between the two lines of cracks has 
been varied. 

Each sample is prepared with the 11 initial cracks forming the array using a cutter 
blade. The experimental set-up consists of a tensile machine driven by a motor (Micro 
Controle UE42) controlled electronically to move step by step (Micro Controle ITL09). 
The paper sheets are mounted on the tensile machine with both ends attached with 
glue tape and rolled twice over rigid bars clamped on jaws. The motor controls the 
displacement of one jaw (400 steps per micrometer) while the other jaw is rigidly fixed 
to a force gage (Hydrotonics-TC). The tensile machine is placed in a box with a humidity 
level stabilized at 5% and at room temperature. The samples are loaded by increasing 
the distance between the jaws on which they are clamped such that the resulting force 
F = acH (e = 50/im is the film thickness and H = 20.95cm the sample height) is 
perpendicular to the initial cracks direction. The loading is enforced at a constant and 
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Fig. 1. Geometry of the crack array initiated in the paper sample and the direction 
of the apphed stress a. 



very slow velocity (between 0.84/im.s~^ and 2/im.s~^) until the final breakdown of the 
sample. More details about the experimental set-up can be found in [7]. Actually, once 
the force threshold F^ of the sample is overcome, we observe a quasi- instantaneous 
rupture of the sample with a quasi-simultaneous growth of all the cracks. There is no 
visual evidence of any sub-critical crack growth process before the final breakdown of 
the sample. During this breakdown, the applied force drops. When it reaches about 
half its maximum value, we stop loading the sample before breaking it in two pieces. It 
is a convenient procedure that allows us afterwards to image properly the crack path 
configurations in each sample. 

3. Post-mortem analysis of crack trajectories 

After each experiment, post-mortem samples are digitized using a high resolution 
scanner. The obtained images are processed in order to get a binary image in which the 
crack paths can be distinguished. In figure [21 we show an example of a digitized sample 
and of the corresponding extracted binary image. This image processing allows to finally 
create for each crack path the profile y{x) that describes its shape as a function of the 
abscissae x along the axis corresponding to the initial crack direction and centered on 
the corresponding initial crack (cf. figure [3]). 



4. Two types of crack path 

For an horizontal spacing d = 0cm between the two crack lines, the initial crack 
array collapses into a single crack line. In figure HI one can see that the crack paths 
show a repulsion phase before an attraction one. This initial crack pattern always leads 
to the so-called type Bq growth behavior. This is the signature of the fact crack tips 
repulse each other. Indeed, the crack tips move away from the initial crack line before 
the neighboring tips overtake each other. Then, the paths curve to join the neighboring 
crack lips tending to form a right angle. Couples of crack paths appear to describe a 
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Fig. 2. Post-mortem images of a fractured sample witli d = 1.5cm : at tlie top, initial 
digitized image and, at the bottom, corresponding binary image with extracted crack 
paths obtained through image processing. 




Fig. 3. On the left, binary image showing the path followed by a crack, and on 
the right, corresponding extracted crack profile y{x). 




Fig. 4. Binary image of a fractured sample with d = 0cm. 



sort of spiral shape as one can see in figure [51 

For < d < (iint (experimentally di^t ~ 1.8±0.2cm), two different kinds of behavior 
have been identified with respect to the crack path shapes among the fifty crack arrays 
that have been fractured. These two behaviors occur statistically for a given sample 
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Fig. 5. Close-up of a binary image of a fractured sample with d = 0cm. 




Fig. 6. Binary image of a fractured sample with d = 0.4cm. 




Fig. 7. Binary image of a fractured sample with d = 5cm. 



geometry. For the first one, that will be referred to as type A, neighboring cracks on 
different lines attract each other from the beginning of the crack growth up to the end 
as it is the case in figure [2] [d = 1.5cm). In contrast, for the type B, the crack paths are 
a mixture of two paths, a repulsive plus attractive path on one end and an attractive 
path tip to tip on the other end of each crack. For instance, in figure [6] {d = 0.4cm), we 
observe a repulsion phase for the right tip of the bottom cracks followed by an attraction 
phase by the neighboring crack on the other line while the right tip of the top cracks 
tends to join the left tip of the bottom cracks in a rather unexpected way. It seems that 
the symmetry breaking associated to the repulsive phase of crack growth on one end of 
the crack allows simultaneously the merging of tips on the other end of the crack. 

Finally, when d > d[^t, the two lines of cracks do not interact (cf. figure [7]). Actually, 
one of the two lines fractures preferentially and we get the type Bq behavior previously 
observed for d = 0cm since the situation is comparable. 

To summarize the results : 

- For d = 0, only type Bq is observed. 

- For < d < dint, types A and B are observed. 

- For dint < d, only type Bq is observed. 

The statistical proportion between types A and B is a function of the spacing d between 
the two initial crack lines. In figure [H we plot for each value of d the percentage P 
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Fig. 8. Percentage of samples for which the crack path are of type B (including the 
case Bq). Experimentally, there is some uncertainty concerning the distance at which 
there is no more interactions between the two crack lines. This is schematically shown 
by the dotted line and lighter color band around d = 1.8cm. 



of samples for which the fracture process follows type B (including the case Bq). This 
percentage decreases regularly starting from 100 % for d = down to zero for d = djnt 
where it jumps back brutally to 100% (this last case corresponds to rupture on a single 
crack line as observed in figure [T]) . 

In the case of a homogeneous material, we would have expected the crack 
trajectories to be reproducible from one experiment to the other. A sheet of paper 
is actually made of a complex network of cellulose fibers. Scanning electron microscopy 
on our samples shows fiber diameters between 4 and 50yU.m with an average of 18/im. 
Cellulose fibers are themselves a bundle of many microfibrils that have a crystalline 
structure [ISj. The statistical selection between two types of crack path shapes for a 
given initial crack geometry is probably triggered by this heterogeneity of the fractured 
material. Indeed, heterogeneity and anisotropy in the initial crack tip toughness and 
shape due to the heterogeneity of the material might induce small perturbations in 
the initiation of the crack growth path. We discover that these perturbations are large 
enough for the crack to explore statistically two very different "meta-stable" crack paths. 

Additionally, the crack patterns observed in this section allow us to confirm the 
fact that, in general, two crack tips repulse each other while a crack tip and a crack lip 
attract each other. 



5. Analysis of attractive cracks in type A samples 

In this section, we study the average behavior followed by type A cracks that 
corresponds to experiments in which the cracks have been attracting each other during 
their whole growth. For each value of the spacing d between the two crack lines, we 
extract the average profile of the cracks. Actually, for each sample, we compute the 
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Fig. 9. (a) Average crack profile {y{x)) for ttiree samples with d = 1.2cni and (b) 
average crack profile {{y(x))) for five values of d, for experiments with a type A 
behavior. 
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Fig. 10. (a) Average crack profile {y{x)) for a sample with d = 1.2cm and its extension 
allowing to define ij and (b) distance ij as a function of the spacing d between initial 
crack lines, for experiments with a type A behavior. Each point corresponds to one 
experiment. 



mean crack profile {y{x)) averaging over tlie profiles y{x) of all the cracks located on the 
same initial crack line. In figure [9](a), we plot for three different samples the averaged 
profiles of the two crack lines separated hj d = 1.2cm. Then, we average the profiles 
{y{x)) over all the type A samples corresponding to the same value of d (cf. figure [9t^b)). 
We get the mean type A profile {{y{x))) for each value of d. 

It is important to notice that, depending on the position x, the number of profiles 
on which the averaging is performed is variable since all the cracks do not have the 
same length in the x direction. This is the reason why the averaged profiles, {y{x)) and 
{{y{x))), are not continuous at all points. For the same value of d, cracks of different 
samples appear to a have very similar mean profiles {y{x)) as one can see in figure [9](a). 
Additionally, for all type A cracks, the dependence of the crack path on the distance 
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Fig. 11. Paths of repelled cracks (type B). (a) Average crack profile {y{x)) for two 
samples with d = Icm, and (b) average crack profile {{y(x))) for four value of d. 



between the initial crack lines d is weak and non-systematic. The main difference is that 
the crack paths extend farther away when the distance d is increased (cf. figure [n](b)). 

Extending intuitively each crack path to join the neighboring crack lip on the other 
crack line, we are able to define an hypothetical junction position ij measured from the 
center of the joined crack (cf. figure [TOT a)). In figure [TOT b). we see that this junction 
distance ij increases rapidly with d. Actually, the rapid increase of this distance between 
d = 1.2cm and d = 1.6cm suggests that there might be a divergence of the junction 
position ij when the distance d approaches (iint(~ 1.8cm). A divergence of ij would 
make sense since when d > rfint the two crack lines do not interact anymore and thus, 
the cracks of one line cannot cross those of the other line. 



6. Analysis of repulsive cracks in type B samples 

In this section, we study in type B (including Bq for d = 0) samples the average 
behavior followed by the repulsive crack paths. In figure [TTT a). we plot the average 
profile of repulsive cracks {y{x)) for two samples with d = 1cm. We note that, for a given 
distance d, the crack path is reproducible from one sample to the other. In figure [TTTb). 
we plot the average crack profile {{y{x))) as a function of d and x. Clearly, the smaller 
d, the larger the repulsion. One way to quantify this is to measure the maximum angle 
of deviation 9max during the repulsion phase. As can be seen in figure fT2l this quantity 
decreases with the distance d. If we extend linearly this curve, we can see that the 
effects of the repulsion disappear when d reaches a value of about 1.9cm. Interestingly, 
we recover a value in the range corresponding to the characteristic distance rfint above 
which the two crack lines do not interact. 

For d = 0, we find a maximum deviation angle Omax — 27°. This value is significantly 
larger than Kachanov's prediction 6max — 14° [32] that was obtained assuming the crack 
goes in the direction for which the strain energy release rate would be maximum. Melin 
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Fig. 12. Maximum angle of repulsion of the crack path during the repulsion phase as 
a function of the distance d for type B experiments. 



[30] has a different approach to explain the deviation of coUinear cracks. She analyzes 
the stability of the crack path to a local perturbation. The crack deviation angle is then 
6 = atan(25y/L) where 6y is the amplitude of the local deviation and L is the crack 
length. To obtain O^ax — 27° would require a perturbation 6y ^ L/4 = 0.25cm, a rather 
unphysical value. A large deviation angle could also be obtained assuming that several 
smaller local deviations of the crack path occur due to the heterogeneous structure of 
paper. However, if that was the case, one would expect the crack to progressively rotate 
instead of turning rather abruptly (see figure [TlTb)). 



7. Discussion and conclusion 

In this article, we have seen that the interactions between two lines of cracks growing 
under a very low strain rate in a heterogeneous material lead to a statistical behavior. 
Indeed, for the same initial sample geometry, we get two types of crack path pattern. The 
type A behavior, that corresponds to a permanent attraction between the cracks during 
their growth, appears to present a universal shape that does not depend significantly 
on the initial crack array geometry. On the contrary, for type B behavior, for which the 
crack paths present either a repulsion phase followed by an attraction or an attraction 
phase tip to tip, the path trajectories appear to be dependent on the distance between 
the two initial crack lines d. Indeed, the repulsion decreases as d increases and finally 
totaly disappears as d tends to a certain distance di^t. This particular value of d is also 
a critical value above which the two lines do not interact anymore. Understanding this 
value is still an open issue. Also, the value of the crack deviation angle observed in our 
experiments during the repulsive phase remains unexplained. In particular, it is rather 
large compared to theoretical predictions. 

The statistics observed for the crack path shape is the signature of the meta- 
stability of two types of crack trajectory for a given geometry. This bi-stability is for the 
moment an unexplained property but it seems to be intrinsic to the chosen crack pattern 
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geometry. It is likely that the trigger for this statistics comes from the complex structure 
of paper sheets that makes this material heterogeneous and leads to anisotropy on the 
initial crack tip shape and toughness. The two types of crack path are two meta-stable 
paths between which each crack has to choose when it initiates its growth. Only at this 
initiation stage, may the heterogeneity of the sample strongly influence the crack in its 
"choice". 

It is also important to point out the collective behavior of the cracks located on the 
same initial crack line. Indeed, inside a given crack line, all the cracks select a specific 
dynamics with quasi-identical crack path shape. Actually, the initiation of the global 
rupture of a sample is probably triggered at one particular crack tip. The choice for the 
crack shape type is made at this particular moment. Then, an unexplained collective 
effect operates to make all the cracks throughout a line behave the same way. 
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